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Motivated by very recent studies of Coulomb drag in grahene-BN-graphene system we develop a 
theory of Coulomb drag for the Fermi liquid regime, for the case when the ratio of spacer thickness 
d to the Fermi wavelength of electrons is arbitrary. The concentration (n) and thickness dependence 
of the drag resistivity is changed from n~''^d~'^ for the thick spacer to n~^ \ \n{nd^) \ for the thin one. 
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Coulomb drag in bilayer semiconductor systems is a 
very interesting phenomenon providing a unique infor- 
mation about many-body effects^""*. Since the role of 
electron-electron interactions in graphene, a novel two- 
dimensional material with extraordinary electronic and 
structural properties^"—, is a controversial issue now 
(for review, see RefHol) study of the Coulomb drag in 
graphene is important, as a way to clarify the situation. 

First theoretical^:-!? and experimentali^ii^ studies of 
the Coulomb drag in graphene have been performed al- 
ready. Theorjiii deals with the case of thick spacer 
{kpd » 1, where kp is the Fermi wave vector of graphene 
and d is the spacer thickness), and the results are in 
a good agreement with the corresponding experimental 
datapi^ (the effects of trigonal warping^? seem to be neg- 
ligible). Very recently, the group of A. Geim and K. 
Novoselov has performed experiments with graphene on 
a substrate of BN, using also BN as a spacer between two 
graphene single layers, and observed the drag for much 
thinner spacersi^. Here I present a theory generalizing 
that of RefflTI for arbitrary values of kpd. 

We start with the general expression for the drag 
conductivitj^i^ based on the lowest-order perturbation 
theory in interlayer Coulomb interaction (we consider 
only the case of identical layers 1 and 2, qualitatively 
all the basic physics remains the same for the case of 
different doping of two layers): 
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where F^; is the nonlinear susceptibility of electrons in 
the layer, Ep is their Fermi energy, T is the temperature, 
Ui2 is the Fourier component of the screened interlayer 
Coulomb interaction which reads, for the case of identical 
layers: 
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Vc{q) and Uc{q) are the Fourier components of bare 
Coulomb interactions within the layer and between the 
layers, respectively, H {q, lo) is the polarization operator 
of the electron gas in graphene. In vacuum, Vc{q) = 
2'Ke'^/q and Uc{q) — Vc{q) exp (—qd) and the expression 
^ coincides with Eq.(A2) from ReflH. 

Since typical frequencies contributing to the integral 
in Eq.([T]) is of the order of ksT/h one can assume, for 
low enough temperatures, that the screening is static and 
replace U12 {q,uj) by U12 (g,0). 

In the ballistic regime (which means that the distance 
between the layers d is much smaller than the mean-free 
path within the layer I) one can calculate for the case of 
grapheneii 
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at a; — > 0, where t — l/v is the mean-free path time and 
V 10^ m/s is the Fermi velocity. 

Assuming that ctd ^ cr where a is the inlayer con- 
ductivity the drag resistivity \s po ~ — ctd/c^. Drude 
formula for the case of graphene is 



cr = -^EpT, 



(4) 



thus, the drag resistivity does not depend on r. 

As a result, the drag resistivity for the case of identical 
graphene layers takes the form 



PD 



2h 
3? 



Ep 



E 



U12 (g,o) 



hv 



(5) 



(2) 



The quantity U12 in Eq.([5]) is given by the expression ^ 
with the static polarization operator for graphene (see, 
e.g., RefsIHIil) 



To find Vi,{q) and udq) one needs to solve the electro- 
static problem taking into account different screening by 
substrate, spacer, and air. Let us assume that the dielec- 
tric medium is three-layer, with the dielectric constant 
distribution 
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The calculations are quite standard (see, e.g., ReflT7|). 
However, for the reader's convenience we present them 
here with some details. 

We have to solve the Poisson equation 



d 
dz 



s{z) 



dip (z) 
dz 



cfe{z)ip{z)^-ATreS{z-ri) (8) 



where if (z) is the electrostatic potential created by the 
point charge e situated at a; = 0, y = 0, z = r/ — > +0. The 
only allowed solution at z < is 



if (z) = v4e«^ 



and at z > d is 
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For rj < z < dit should be tried in the most general form: 

(z) = ae«^ 4- /3e"9^ (11) 

From continuity of the potential and the normal compo- 
nent of electric induction, £)„ — — £^ at the boundaries 
z — Q and z — d we find, taking into account Eqs. Q 
and (fTOl): 
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where prime means the derivative with respect to z. The 
potential (p (z) is continuous at z = 77 {ip{—0) = ip{+0)) 
but its derivative has a jump. Due to Eqs. ([SJ, ©, 
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At last, we can find from Eqs. (ITil) and (IT5|) the constants 
a and /3. The final answer for Vc{q) = eip{z = 0) and 
Uc{q) = eip{z ~ d) reads: 
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For simplicity, we will consider further only the case 
£1 = £2 (which, actually, takes place in the experimental 
situatioi>i^ where BN is used both as a substrate and as 
a spacer). In this case, the expression (ITSt is simplified 
dramatically: 



Uc (q) = Vc (q) exp {-qd) , 
Vc{qj = 
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Substituting Eqs.®, ((15]), ^ into Eq.® we have: 
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where 
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is the effective "fine structure" constant (for the combi- 
nation of BN with £2 ~ 4 and air with £3 = 1 we have 
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FIG. 1: Drag resistivity (|18[) as a function of charge carrier 
concentration, for £3 = 1, £2 = 4, T = 120 K, d — 3 nm 
(dashed line) and d — 4 nm (solid line) 



a w 0.87) and the function F is represented as 
00 
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and Eas.p8|. (l20l) give the known resultii 
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where k = Aakp is the inverse screening radius. 

In the opposite limit y ^ 1 typical values of a; w 1/y ^ 
1 and one can assume (^{x) ~ irx/A which gives 
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The behavior of drag resistivity as a function of charge 
carrier concentration for kpd of the order of one is shown 
in Fig. 1. This result seems to be in a qualitative agree- 
ment with the experimental dataM, at least, it gives the 
correct order of magnitude for the drag resistivity. At the 
same time, for small enough kpd and a « 1 the interlayer 
Coulomb interaction is in general not small, and it is not 
clear whether the lowest-order perturbation theory used 
here will be also quantitatively accurate or taking into 
account next-order contributions will be necessary. The 
issue requires further studies, both experimental and the- 
oretical. 

Recently two more works on the subject appeare d^^i^^ . 
The results of this paper and the other two papers 
concerning concentration dependence of the drag re- 
sistivity are in an agreement, namely, in RefllSl the 
same analytical concentration dependence as here, esc 
n~^\\n{nd'^)\, was obtained for the case of thin spacer 
whereas in ReflU the numerical data were fitted by 
Pd oc n^" with a of the order of one. 
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